Applied Energy 92 (2012) 395-405 



ELSEVIER 


Contents lists available at SciVerse ScienceDirect 

Applied Energy 

journal homepage: www.elsevier.com/locate/apenergy 



Optimization of geometrical parameters for Stirling engines based 
on theoretical analysis 

Chin-Hsiang Cheng*, Hang-Suin Yang 

Institute of Aeronautics and Astronautics, National Cheng Kung University, No. 1, Ta-Shieh Road, Tainan 70101, Taiwan, ROC 


ARTICLE INFO 


ABSTRACT 


Article history: 

Received 13 March 2011 
Received in revised form 23 August 2011 
Accepted 14 November 2011 
Available online 13 December 2011 


Keywords: 

Stirling engine 
Optimization 
Theoretical analysis 
Geometrical parameters 


This study is aimed at theoretical analysis of the effects of the geometrical parameters on the shaft work 
of the Stirling engines. The optimal combination of the phase angle and the swept volume ratio, that leads 
to maximization of the shaft work of the engine, is obtained under different specified conditions. Effects 
of the effectiveness of mechanism, the dead volume ratio, and the temperature ratio on the maximum 
shaft work of the engine as well as the optimal combination of the phase angle and the swept volume 
ration are evaluated. Theoretical analysis of the performance of three types of Stirling engines, at-, p-, 
and y-type, has also been carried out, and a comparison in relative performance among these three types 
of engines is attempted. Results show that for the particular cases considered in this study, the 
p-type Stirling engine produces highest shaft work and the y-type engine the lowest. In general, the y- 
type engine must be very mechanism effective so as to deliver sufficient shaft work. However, among 
the three types of engines, the y-type engine is most capable of operating with low temperature differ¬ 
ence. On the contrary, the a-type engine is particularly not suitable for the applications with low temper¬ 
ature difference since its dimensionless shaft work is found to gradually vanish as the temperature ratio 
is increased. 

© 2011 Elsevier Ltd. All rights reserved. 


1. Introduction 

Stirling engines are referred to as external combustion engines, 
and hence, they can be operated with a variety of heat sources and 
have been applied for energy conversion in a number of engineer¬ 
ing applications [1-3]. The ideal regenerative Stirling cycle consists 
of two isothermal processes and two isochoric processes, as shown 
in Fig. 1. In theory, the thermal efficiency of an ideal regenerative 
Stirling cycle is equal to the Carnot cycle efficiency under the same 
thermal reservoirs temperatures. In practices, the Stirling engines 
are classified into three types in terms of configuration [4j. The 
schematics of the three types of Stirling engines and their varying 
displacements of the pistons/displacers are plotted in Fig. 2. In an 
oc-type engine, the total volume of the working fluid in the engine 
is determined by the relative displacements of its two pistons, 
whereas in a p- or y-type engine, the total volume of the working 
fluid in the engine is simply determined by the displacement of 
piston and is not influenced by the position of displacer. 

According to the literature survey, the first successful analytical 
model for the Stirling engine was presented by Schmidt [5] in 
1871. However, based on the Schmidt model, a group of research¬ 
ers, for example, Finkelstein [6], Walker [7], and Kirkley [8], inves¬ 
tigated the optimal phase angle and swept volume ratio for the 
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oc-, p-, and y-type Stirling engines and refrigerators. It appears in 
these previous studies that among these three types of engines, 
the P-type Stirling engine produces highest indicated work 
whereas the y-type the lowest. 

Lately, the non-isothermal condition was introduced to the anal¬ 
ysis of Stirling cycle by Finkelstein [9]. In recent years, theoretical 
analysis was extended to the irreversible cycle analysis [10-14], 
and some numerical models for Stirling engine were also developed. 
Martaj et al. [14] presented a thermodynamic analysis of a low tem¬ 
perature Stirling engine at steady state operation, and energy, entro¬ 
py and exergy balances were presented at each main element of the 
engine. Recently, Timoumi et al. [ 15,16] and TIili et al. [17] optimized 
the regenerator and swept volume of the Stirling engine by using 
numerical model. Cheng and Yu [18] proposed a thermodynamic 
model for the P-type Stirling engine which taking into account the 
effectiveness of the regenerative channel as well as the thermal 
resistances of the heating and cooling heads. Lately, the thermody¬ 
namic model was then incorporated with a dynamic model to per¬ 
form a dynamic simulation model for the P-type Stirling engines 
with cam-drive mechanisms by the same group of authors [19]. 
More recently, Yang et al. [20] applied thermoacoustic theory to ana¬ 
lyze a free-piston Stirling engine. 

Note that the regenerator effectiveness has been recognized as 
an important factor to the performance of the Stirling engines. 
The regenerator effectiveness is dependent on porosity, permeabil¬ 
ity and material of porous matrix. Temperature distribution in the 

















396 


C.-H. Cheng, H.-S. Yang/Applied Energy 92 (2012) 395-405 


Nomenclature 



E 

effectiveness of mechanism 

W_ 

forced work (kj) 

m 

mass of working fluid (kg) 

W 

dimensionless indicated work 

p 

pressure of work space (kPa) 

W s 

dimensionless shaft work 

Pac 

average cyclic pressure, = ^p m31l p min , (kPa) 


dimensionless forced work 

Pb 

back pressure (kPa) 



P mean 

average pressure, =(p max + p min )/2 (kPa) 

Greek 

Symbols 

Qin 

dimensionless heat input 

a 

phase angle (°) 

R 

gas constant (kj/(kgl<)) 

n 

thermal efficiency 

r 

compression ratio, =V max /Vmin 

K 

swept volume ratio, =V cs /V es 

To 

temperature of compression space (K) 

T 

temperature ratio, =T C /T e 

T d 

temperature of dead space, =(T e + T c )/2 (I<) 

<t> 

crank angle (°) 

Te 

temperature of expansion space (I<) 

X 

dead volume ratio, =V d /V es 

V 

volume variation of total working space (m 3 ) 



Vc 

volume variation of compression space (m 3 ) 

Subscripts 

Vcs 

piston swept volume (m 3 ) 

ac 

cyclic average 

V d 

volume of dead space (m 3 ) 

c 

compression space 

V e 

volume variation of expansion space (m 3 ) 

d 

dead space 

Ves 

displacer swept volume (m 3 ) 

e 

expansion space 

V r 

volume of regenerator (m 3 ) 

max 

maximum 

Vj 

total swept volume, =17 CS + V es (m 3 ) 

mean 

average 

w 

indicated work (kj) 

min 

minimum 

w s 

shaft work (kj) 

r 

regenerator 



Fig. 1 . Ideal regenerative Stirling cycle. 

regenerator can be measured by experiment to determine the 
regenerator effectiveness. Furthermore, the magnitude of the 
regenerator effectiveness is also dependent on the operating speed. 
Actually, they decrease with an increasing operating speed. This is 
because the time for heat exchange between the gas and the por¬ 
ous matrix of the regenerator is decreased by increasing the oper¬ 
ating speed. 

In addition to the above thermodynamic models, a three- 
dimensional computational fluid dynamics (CFDs) approach may 
also be adopted in the analysis of the Stirling engines. The CFD sim¬ 
ulation of a solar-powered Stirling engine was performed by Mah- 
kamov [21]. Numerical predictions of the flow, pressure and 
thermal fields in a reciprocating piston-cylinder assembly of a Stir¬ 
ling engine were carried out by means of the CFD methods by 
Cheng and Hung [22]. It has been recognized that the CFD methods 
give a certain amount of detailed information that are not easily 
obtained by the experiments. However, the three-dimensional 
CFD computation is in general rather time consuming and requires 
much more computation efforts than the thermodynamic models. 


On the other hand, experiments have also been conducted by a 
group of researchers. To name a few, Batmaz and Ustun [23] de¬ 
signed an oc-type Stirling engine and measured the power at differ¬ 
ent heating temperature and charged pressure. Cinr et al. [24] and 
Karabulut et al. [25] performed experimental study on the devel¬ 
opment of a [3-type Stirling engine working at different charged 
pressure and heating temperature. Kongtragool and Wongwises 
[26] conducted experiments for the low temperature differential 
y-type Stirling engines. In these experimental studies, the engine 
torque, the shaft work, and the thermal efficiency at different en¬ 
gine speeds were measured. 

Senft [27,28] considered the effect of the mechanical loss and 
proposed a concept of the mechanism effectiveness. The mecha¬ 
nism effectiveness is well known as a ratio measure of how well 
a mechanism utilizes work if it is regarded as a work transmitter. 
The value of mechanism effectiveness can be estimated by using 
the dynamic model similar to that proposed by Cheng and Yu 
[19]. As stated in Refs. [27,28], the typical value of the mechanism 
effectiveness is within the range between 0.7 and 0.9. In the con¬ 
cept, the shaft work of engine can be determined in terms of the 
indicated work and the forced work. In addition, the indicated 
work is calculated simply by integrating the area inside the cycle 
in a p-V diagram and the forced work is determined by calculating 
the areas of the shaded zones shown in Fig. 3, where p b is the back 
pressure exerted on the opposite side of the piston. The back pres¬ 
sure could usually be the atmospheric pressure or an elevated 
pressure in case a pressurized space is placed on the opposite side 
of the piston as a back zone. 

The fundamental problems that the designers of the Stirling en¬ 
gines have to face are the time consumption in the design process 
and the cost of precision manufacturing. In most of the cases, the 
determination of the geometrical parameters of the engines is based 
on the experience of the designers or on insufficient costly experi¬ 
mental information. Therefore, the improvement of the engine can 
only be carried out in a relatively high-cost, time-consuming 
trial-and-error process. In this regard, modeling the engines is ex¬ 
pected to be one of the alternative solutions to these issues. Further¬ 
more, it is important to mention that the shaft work is exactly the 
output work delivered by the engine. However, in most of these 
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Fig. 2. a-, p-, and y-type Stirling engines and their varying displacements of pistons and displacers. 



Fig. 3. Indicated work and forced work [27,28], 

existing models [5-22] it is the indicated work that was maximized, 
not the shaft work. Therefore, the performance of the engine might 


have not been fairly assessed. Under these circumstances, the pres¬ 
ent study is aimed at optimization of the geometrical parameters for 
the Stirling engines by introducing the concept of the mechanism 
effectiveness into the theoretical model. The optimal combination 
of the phase angle and the swept volume ratio, that leads to maximi¬ 
zation of the shaft work delivered by the engine, is then investigated. 
Meanwhile, a comparison in the relative performance of the a-, p-, 
and y-type engines is attempted. 

2. Theoretical analysis 

Here, the reversible model is adopted. A simple idealized model 
enables researchers to study the effects of major parameters which 
dominate the cycle without getting bogged down in the details. 
The present study is somewhat idealized, but they still retain the 
general characteristics of the real cycles they represent. 

In this study, a group of dominant parameters are taken into ac¬ 
count. They are: mechanism effectiveness (£), phase angle (a), 
swept volume ratio (k), dead volume ratio (x), and temperature ra¬ 
tio (t). The first four parameters are referred to as the geometrical 
variables, and the temperature ratio is an operating parameter. 
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During the operation of the Stirling engine, the volume varia¬ 
tion of the expansion space should move ahead of the compression 
space by a leading time. Hence, there exists a phase angle (a) be¬ 
tween the movements of expansion and compression pistons (a 
type) or between displacer and piston (p and y types). 

In addition, the swept volume ratio, the dead volume ratio and 
the temperature ratio are defined as 

K = V a /V B , x = V d /V es , and r = T C /T e (1) 


respectively, where V cs , V es , and V d denote the piston swept volume, 
the displacer swept volume, and the dead volume. The dead volume 
(V d ) of the engine is determined by V d = V em j n + V cmirl + V r , where 
V emin and Vanin are the clearance volumes of the expansion and com¬ 
pression spaces, respectively, and V r is the volume of regenerator 
channels. 

Dimensionless shaft work, dimensionless indicated work and 
dimensionless forced work are defined, respectively, as 


W s = 


Ws 

mRT e ' 


W = 


W 

mRT e ' 


W_ 


W_ 

mRT e 


( 2 ) 


The dimensionless parameters are based on hot temperature and 
are different from those using average pressure and volume, and 
hence they are more accurate. When the dimensionless shaft work 
is fixed, more mass of the working fluid or higher hot temperature 
means higher shaft work. Note that the dimensionless shaft work 
(W s ) is the objective function to be maximized in the present study. 

Theoretical analyses of different types of Stirling engines are de¬ 
scribed individually in the following: 

In an a-type Stirling engine, the volumes of the expansion and 
the compression spaces can be expressed as 

Ve(<P) = V e ^+lv es (\+ sin(4» + oc)) (3a) 


Vc(<l>) = V c ^+^V es K( 1+sin^) (3b) 

For the p-type or the y-type Stirling engine, the relation of the 
volume variation of the expansion space is similar to that for the 
a-type engine, which is dependent on the displacement of the dis¬ 
placer only. However, the volume of the compression space is 
dependent on the displacements of both displacer and piston. That 
is, for the p-type engine: 

VeW = \4 mln +iv«,(l+sin(* + a)) (4a) 


=x + ^( 2 + k) +if<sin0 for y type 


(7c) 


The compression ratio of the engine is determined in terms of 
the maximum and the minimum values of the total volume as 


r _ ^max _ 1 | 2\/l+2K-COSa+K- 2 

^min 2X+1 + K- yi+2KCOSa+K 2 


r _ ^max _ i _2k_ 

Vmin 2^+1 -K + Vl -2KCOSOC+K 2 


for a type 
for p type 


(8a) 

(8b) 


r = = 1 + TTv for 7 type ( 8c ) 

According to the information provided in Refs. [13,18,19,28], 
the pressure drop between the working spaces in the Stirling en¬ 
gine is negligible. Therefore, the pressure is assumed to be uniform 
throughout the engine and can be calculated in terms of the vol¬ 
umes and the temperatures of the working spaces as 

mR 

P ~ v, i Vc_ , Yi 

T e + T c + T d 

where the temperature of dead space is approximated by 
T d = (T c + T e )/2. The pressure can be further expressed in the fol¬ 
lowing form: 

mRT e 


r V , 
where 

W(P>,(x,K,r,x) 


W{ 4 >, a,K,x,x) 


2x K + t Ksiiuj) sin(<jS + a) 


L(1 + T) 


2r 


2t 


( 10 ) 


for a type (11a) 


¥(<!>,01,K,T,X) 


2x + T + yi -2fv'cosa + fc 2 + Ksin0 + (T-l)sin(<i>+a) 


(1+T) 


2t 


2t 


2t 


for/ltype (lib) 


W(<j>,a,K,z,x) 

2x 1+t + k Ksin (</>) (t- l)sin(</> + a)l" 


(1+T) 


2t 


2t 


2t 


J 


for y type 


(lie) 


Combining Eqs. (2), (6), (7a), (10) and (11) leads to the dimen¬ 
sionless indicated work as 


V c {(j>) = W min +^es(\/l - 2fccosa + k 2 + Ksin0- sin(<£ + a)) 

(4b) 

and for the y-type engine: 

V e (4) = V e ^+} i V es (\+sm(4> + a)) (5a) 

VcW = W mm + ^ V K (1 + k + K sin - sin(0 + a)) (5b) 

Then, the total volume ( V) can be calculated with 


W = 


fpdV 


l 


2n d(f> . 
'V ,+-d</> 


mRTe Jo * d<t> 

And the dimensionless heat input can be calculated by 
fpdV e 


Qin = 


mRT e 


1 r 

2 J 0 


W cos(4> + a)d 4> 


( 12 ) 


(13) 


Integrating the terms on the right-hand side of Eqs. (12) and 
(13) by using residual theorem, one obtains the dimensionless 
indicated work and dimensionless heat input of the three types 
of Stirling engine as 


V(4) = V e {<t>) + VcW + V r = Ves<P(<t>) 

where 


( 6 ) 


<P{4>) = Z + ^(l + K) + ^(sin(<£ + a) + zesin0) for a type (7a) 


W = 



27tiK(l - t) sin a 



n 27tTKsina 

(p- 


(^1 + v 1 -2kcosoc+k 2 ) +-/csin</> for/?type 

(7b) a i + b 2 


Jp 2 -(a 2 +b 2 ) J 


(14a) 


(14b) 
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Fig. 4. Effects of swept volume ratio on dimensionless works, at E = 0.8, z = 0.3, and 
Z-0.1. 

where for ot type: 
a = isina 

b = K+ Tcosa (15a) 

P = 4t%/( 1 +T) + K + T 

for (3 type: 
a = -(1 - t) sin a 

b= k- (1 - r)cosa (15b) 

p = 4t^/(1 + t) + t + v 1 -2fccosa + K 2 


Fig. 5. Effects of phase angle on dimensionless works, at £ = 0.8, x = 0.3, and x = 0.1. 

for y type: 

a = -(1 - t) sin(a) 

b = k — (1 - t) cos(a) (15c) 

p = 4xyJ(\ +T) + 1+ T + K 

Note that the magnitudes of the forced work and the indicated 
work are dependent on the back pressure (p b ). As discussed by Sen- 
ft [29], the optimal value of the back pressure is the average cyclic 
pressure, which is 


Pb \/PminPmax 


( 16 ) 
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(a) Effects of dead volume ratio, at 
E=0.8, a= 80°, k= 0.5, and r=0.3. 



(b) Effects of temperature ratio, at 
E=0.8, a=80°, k=0.5, and x=0.1. 


Fig. 6. Effects of dead volume ratio and temperature ratio on dimensionless works 
for a-type engine. 



(a) (3 type with E=0.8, a=30°, k= 1.2, and %=0.1. 



(b) y type with E=0.8, a=100°, k=1.0, and y=0.1. 

Fig. 7. Effects of temperature ratio on dimensionless works for )i- and 7-type 
engines. 

2.1. a Type 


where p min and p max are the minimum and maximum values of the 
pressure in a cycle. Introducing Eqs. (10) and (11) into the above 
equation leads to 


Both the two pistons in the expansion and the compression 
spaces of the a. type Stirling engine keep contact with the back 
pressure; therefore, one needs to consider the forced works done 
by the two pistons, which are 


Pb = —p—^ V'^min’f'max (17) 

v es 

where 

^y min 'F max = 2T ( 18 ) 

\/P 2 - ( q2 + b 2 ) 

One can calculate the dimensionless shaft work for the engines 
as described by [28,29], which can be written in dimensionless 
form as 


l~ E ) w - 0 9 ) 

where 

W_ = W e _ + W c _ (20) 

and the dimensionless forced works done by the pistons in the 
expansion and the compression spaces for different types of engines 
can be determined as follows: 



We_= j>[(p-p b ) dV e ] = j [( S' - yf w min T 7 max) COS {<j> + «)d</>] 

(21a) 

= j [(p - Pb )dV c ]~ = TEl j [( W - V'l'm.n'/'max) KCOSc/,df/>] 

(21b) 


where the superscripted negative sign means 

[„] 1 = / ° if n >° 

U \n if n < 0 

In dimensionless forms, the above equations become 

TT7 1 


W e _ = ^ f |(7'- 


min T 7 max) COS(<)> + CC)d(j) 


(21c) 


(22a) 



[(T 7 - '/T'minT'max) K COS <£d</>j 


(22b) 


2.2. p Type 

One of the major differences between p and ot types is that the 
p-type engine has only one piston that keeps contact with the back 






















C.-H. Cheng, H.-S. Yang/Applied Energy 92 (2012) 395-405 


401 



Fig. 8. Dimensionless shaft work as a function of phase angle and swept volume 
ratio for oc-type engine at E = 0.8, x = 0.3, and % = 0 . 1 . 

pressure. Therefore, only a forced work done by the piston needs to 
be considered. 

W_ =W c _=^j [(S'- \/ f'jnin COS (23) 

2.3. y Type 

The calculation of forced work for y type is similar to that for p 
type since the y-type engine also has only one piston in the com¬ 
pression space that keeps contact with the back pressure. Thus, 
the dimensionless forced work of y type Stirling work can be calcu¬ 
lated with Eq. (23). 

3. Results and discussion 

A parametric study is performed based on the theoretical anal¬ 
ysis so as to evaluate the effects of the dominant parameters. Fig. 4 
shows the effects of the swept volume ratio on the dimensionless 
works, at E = 0.8, t = 0.3, and % = 0.1, for all types of engines. Plot¬ 
ted in Fig. 4a are the results of W, W s , and W_ for a-type Stirling 
engine with the phase angle a = 80°, Fig. 4b for p type with 
a = 30°, and Fig. 4c for y type with a = 100°. It is clearly seen that 
for all types of engines the dimensionless shaft work first increases 
with the swept volume ratio in the 1ow-k regime. When the swept 
volume ratio is elevated, the shaft work reaches a maximum at a 
critical value of k. Over the critical swept volume ratio, the shaft 
work gradually descends in the high-K regime. The curves of 
dimensionless indicated work may also reach a maximum value; 
however, the peaks of the curves of the shaft work and the indi¬ 
cated work take place at different swept volume ratios. This im¬ 
plies that the optimal swept volume ratio obtained in terms of 
the indicated work is different from that in terms of the shaft work. 


Note that the peak shaft work and its accompanying critical swept 
volume ratio are dependent on the specified values of a, £, t, and x- 

Fig. 5 conveys the effects of the phase angle on the dimension¬ 
less works, at E = 0.8, t = 0.3, and /_ = 0.1. Fig. 5a shows the results 
for a-type Stirling engine with the swept volume ratio k = 0.5, 
Fig. 5b for p type with k = 1.2, and Fig. 5c for y type with k = 1.0. 
Again, a peak value of the dimensionless shaft work is observed 
for each of the a-, p-, and y-type engines. In the three plots of this 
figure, it appears that the critical phase angle is approximately 80° 
for a type, 30° for p type, and 100° for y type. These obtained opti¬ 
mal phase angles agree closely with the results provided by Kirkley 
[8] and Senft [28]. It is also noticed that for the particular cases 
considered in this study, the p-type Stirling engine produces high¬ 
est shaft work and the y-type engine the lowest. 

Shown in Fig. 6 are the respective effects of the dead volume ra¬ 
tio and the temperature ratio on the dimensionless works of the a- 
type engine. Fig. 6a displays the influence of dead volume ratio at 
E = 0.8, a = 80°, k = 0.5, and t = 0.3, and Fig. 6b shows the influence 
of temperature ratio, at £ = 0.8, a = 80°, k = 0.5, and / = 0.1. It is ex¬ 
pected that as the dead volume in the engine is increased, the out¬ 
put work of the engine is reduced. Fig. 6a clearly reflects this 
expectation. It is observed that as the dead volume ratio is in¬ 
creased, the dimensionless indicated work and shaft works are 
both monotonically decreased. 

In addition, as the temperature ratio is elevated, the tempera¬ 
ture difference between the expansion and the compression spaces 
is actually reduced and as a result, the indicated work and the shaft 
work will be decreased. In Fig. 6b, it is noted that as the tempera¬ 
ture ratio approaches 0.63, both the curves of W and experi¬ 
ence a sudden drop to zero, and meanwhile the dimensionless 
shaft work is found to gradually vanish before the temperature ra¬ 
tio approaches a limit of 0.63. This indicates that the a-type engine 
is inherently not suitable for the applications with low tempera¬ 
ture difference. Note that the limiting value of the temperature ra¬ 
tio for the a-type engine is dependent on the £ value. 

Effects of temperature ratio on dimensionless works for p- and 
y-type engines are displayed in Fig. 7. The cases considered in 
Fig. 7a are the P-type engines with £ = 0.8, a = 30°, tc = 1.2, and 
X = 0.1, and in Fig. 7b the y-type engines with £=0.8, a = 100°, 
k= 1.0, and x = 0.1. it is observed that the dimensionless works 
of the two types of engines are monotonically decreased with an 
increasing temperature ratio. Furthermore, as stated earlier, the 
a-type engine is not able to deliver shaft work at a temperature ra¬ 
tio elevated to 0.63. Nevertheless, in Fig. 7a and b, it is observed 
that for the p-type engine when the temperature ratio is elevated 
to 0.63, the shaft work still stays at approximately 33% of the work 
at t = 0.3, and the y-type engine stays at 54%. This indicates that 
among the three types of engines, the y-type engine is most capa¬ 
ble of operating with low temperature difference even though its 
shaft work output is relatively low. 

It is noticed that the temperature ratio t represents the ratio of 
the temperature in compression space to that in expression space. 
Actually, the operating temperature ratio is dependent on the ther¬ 
mal resistance of heating and cooling spaces, the effectiveness of 
regenerator, the working fluid mass, and the engine speed [18]. 

In accordance with the parametric study discussed above, one 
may find that the shaft work of the engine is elevated when the 
dead volume ratio or the temperature ratio is decreased. How¬ 
ever, the relation between the shaft work and the dead volume 
ratio or the temperature ratio is monotonical and hence, there is 
no need to optimize the dead volume ratio or the temperature 
ratio. On the other hand, the peak values of the shaft work 
can be observed as the phase angle or the swept volume ratio 
is varying. Therefore, optimization task of this study is focused 
on the two parameters. Fig. 8 conveys the dimensionless shaft 
work as a function of phase angle and swept volume ratio for 
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(a) £=0.7 


(b) £=0.8 


Fig. 9. Dependence of optimal W s , a, and k on specified E, t, and / for 7-type engine. 


the a-type engine at £ = 0.8, r = 0.3, and x = 0.1. It is found that 
the dimensionless shaft work reaches the maximum of 0.343 at 
a = 80.5° and k = 0.504. However, it is important to mention that 
the optimal combination of the phase angle and the swept vol¬ 
ume ratio and the maximum value of the dimensionless shaft 


work are dependent on the specified conditions specified in 
terms of mechanism effectiveness, temperature ratio, and dead 
volume ratio. Once the specified condition is changed, the opti¬ 
mal geometrical parameters combination that results in the 
maximum shaft work is altered. 
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(a) £=0.7 

Fig. 10. Dependence of optimal W s , a, and 





(b) £=0.8 

on specified E, z, and / for p-type engine. 


Fig. 9 conveys the dependence of the optimal values of the 
phase angle and the swept volume ratio and the maximum of 
the dimensionless shaft work on E, t, and % for the a-type engine. 
It is found that as the mechanism effectiveness is increased from 
0.7 to 0.8, the maximum value of the dimensionless shaft work is 
nearly doubled. The dimensionless shaft work maximum is 


decreased by increasing the value of t or x■ Meanwhile, the opti¬ 
mal phase angle is elevated from 80° to 180° as the temperature 
ratio is varied from 0.3 to 0.6. In general, the optimal phase angle 
is decreased by reducing the value of t or x or by improving the 
mechanism effectiveness (£). The influence of x on the optimal 
phase angle is more remarkable when £ is higher. In addition, 
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(a) £=0.8 


Fig. 11. Dependence of optimal W s , a, and 


the effects of the temperature ratio on the optimal swept volume 
ratio are more involved. At a lower dead volume ratio, say less than 
0.3, the optimal swept volume ratio increases with the tempera¬ 
ture ratio. However, as the dead volume ratio exceeds 0.5, the opti¬ 
mal swept volume ratio may be decreased by increasing the 
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(b) £=0.9 

on specified E, t, and / for y-Lype engine. 


temperature ratio. It is noted that when the mechanism effective¬ 
ness is improved, the optimal swept volume ratio is also increased. 

Dependence of the optimal values of W s , a, and k on specified E, 
t, and x for the |3-type engine is illustrated in Fig. 10. Basically, the 
trend of the curves plotted in this figure is similar to those in Fig. 9. 
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The optimal shaft work delivered by the p-type engine is generally 
higher than that by the ot-type engine. However, the optimal swept 
volume ratio is decreased by increasing the temperature ratio for 
all values of Again, it is found that when the mechanism effec¬ 
tiveness is improved from 0.7 to 0.8, the performance of the engine 
is greatly increased. 

Fig. 11 shows the optimal W s , a, and k at different E, t, and % f° r 
the y-type engine. In this figure, the mechanism effectiveness is as¬ 
signed to be 0.8 and 0.9; however, it is noticed the optimal shaft 
work of the y-type engine is still much less than those of the al¬ 
and the P-type engines shown in Figs. 9 and 10. This seemingly 
means that the y-type engine must be very mechanism effective 
so as to deliver sufficient shaft work. 

4. Concluding remarks 

Theoretical analysis has been performed for the effects of the 
geometrical parameters on the shaft work of the Stirling engines. 
The optimal combination of the phase angle and the swept volume 
ratio toward maximization of the shaft work of the engine is inves¬ 
tigated under different specified conditions. Effects of the effective¬ 
ness of mechanism, the dead volume ratio, and the temperature 
ratio on the maximum shaft work of the engine as well as the opti¬ 
mal combination of the phase angle and the swept volume ratio are 
evaluated. Theoretical analysis of the performance of three types of 
Stirling engines, oc-, p-, and y-type, has also been carried out, and a 
comparison in relative performance among these three types of en¬ 
gines is completed. 

In accordance with the obtained results of this study, it is found 
that the P-type Stirling engine produces highest shaft work and the 
y-type engine the lowest. In general, the y-type engine must be 
very mechanism effective so as to deliver sufficient shaft work. 
However, among the three types of engines, the y-type engine is 
most capable of operating with low temperature difference. On 
the other hand, the oc-type engine is particularly not suitable for 
the applications with low temperature difference since its dimen¬ 
sionless shaft work is found to vanish at a higher temperature 
ratio. 
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